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Abstract 

A number of relations between the Kaplan- Yorke dimension, phase space contraction, transport 
coefficients and the maximal Lyapunov exponents are given for dissipative thermostatted systems, 
subject to a small external field in a nonequilibrium stationary state. A condition for the extensivity 
of phase space dimension reduction is given. A new expression for the transport coefficients in 
terms of the Kaplan- Yorke dimension is derived. Alternatively, the Kaplan- Yorke dimension for a 
dissipative macroscopic system can be expressed in terms of the transport coefficients of the system. 
The agreement with computer simulations for an atomic fluid at small shear rates is very good. 
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1. INTRODUCTION 

The Kaplan- Yorke (KY) or Lyapunov dimension was introducedS as a conjecture relating the HausdorfF 
(H) dimension and the Lyapunov exponents of the invariant measure of a given dynamical system.]^ This 
allows a computation of the H-dimension of the attractor of a dissipative dynamical system in^erms 
of its Lyapunov exponents. Its validity has been proven for two dimensional dynamical systemsIS and 
for a rather large class of stochastic systems. B In general one only knows that the KY-dimension is 
an upper bound for the H-dimension.a Although it is not hard to construct rather artificial counter 
examples, it is generally believed that the conjecture holds for "generic" dynamical systems. 



In this paper we will discuss a relation (see Eq. (4.10) below) between the KY-dimension in large 
thermostatted systems and their physical properties such as the transport coefficients in a nonequilibrium 
stationary state. Such a relation allows us to estimate the KY-dimension from a measurement of the 
transport coefficient. In doing so we obtain a new relation between a dynamical quantity (the KY- 
dimension) and a physical quantity (the transport coefficient). 

A difficulty in doing this is that while the dynamical quantities are usually defined for any finite 
number of particles TV, the physical quantities usually refer to systems of very large N, so that one 
can meaningfully define intensive quantities, which only depend on intrinsic parameters like the number 
density n — N/V rather than on N and V separately, where V is the volume of the system. We would 
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like to state this by saying that strictly speaking, a thermodynamic limit has to be taken, i.e. iV — s- cx), 
V —> oo with the number density {N/V — > n) and other intensive quantities, such as in particular 
the shear rate, considered in this paper, held constant. This is straightforward if the linear transport 
coefficients are required and the limit of the external field — > is taken before the limit — » oo, as 
in linear response theory. However, for finite fields, changes in the behavior of the system can occur, 
when the thermodynamic limit is approached - as e.g. the onset of turbulence in a sheared system. B 
In this paper we are interested only in the behavior of systems before such a transition takes place, like 
the laminar flow of a sheared fluid considered in section 5. We think, nevertheless, that our results can 
be usefully formulated for large systems using expr essi ons like "for sufficiently large N" without taking 
the mathematical limit (see comment after Eq. (2^) for a more precise discussion). Although this 
expression is not mathematically well defined, we think that its meaning will be clear in any practical 
application (see note 7 on page 14 for an attempt to clarify this point). Moreover, because errors in 
intensive quantities are typically 0{N~^) where N ~ 10^^, the approach is physically reasonable. 

The above mentioned connection between a dynamical and thermodynamic treatment requires the 
usually discrete Lyapunov spectrum to be effectively replaced by a continuous intensive spectrum and 
an intrinsic version of the KY-dimension to be introduced. In section 4 we show how this can be 
implemented, after having introduced the basic equations which connect the dynamical and physical 
quantities in section 2, and deriving a new exact relation for the linear transport coefficients in section 
3. 



2. BASIC RELATIONS AND DEFINITIONS 

As has been shown before, % there is a direct relationship between the sum of the non-zero Lyapunov 
exponents Xi with Xi > A^+i, 1 < i < 2dN — f (where N is the number of particles, d the Cartesian 
dimension and / the number of zero Lyapunov exponents) and the phase space contraction rate in a 
thermostatted system, subject to an external force F,,, in a non equilibrium stationary state, of the 
form: 

2dN-f 

J2 A,,jv(i^e) = Ajv(Fe). (2.1) 

4=1 

Here the subscript N indicates the A^-dependence of the various quantities in Eq. (|2.l|), 2dN — / is 
the efl'ective number of degress of freedom in phase space of the system and A^iFe) = ^ • F is the 
phase space contraction rate, where F stands for the collection of the coordinates and momenta of the 
A'' particles and F for its time derivative. For macroscopic systems, i.e. systems with very large N, 
can use the equality of the dynamical phase space contraction and the physical entropy production 



and obtain from Eq. (2.1) 



1 Y^^^ <^N[Fe) _ JN{Fe)Fe _ LN{Fe)F^ 

N ^ ^''"'^^"^ - - nksT (^-^^ 

i—l 

Here CTjv(fe)j "^Ar(-Fe) and Li\[(Ff,) are the entropy production rate per unit volume, the dissipative flux 
and the transport coefficient respectively, induced in the system in the stationary state by the external 
force Fe, where a nonlinear constitutive relation J^iFe) = —Lpf{Fe)Fe has been used. The subscript N 
indicates the A^-dependence for finite systems. The kinetic temperature T is determined by the relation: 



,^_-7^E^-^-^^ (2-3) 



1=1 

where m is the particle mass, {p;, i = 1, N} are the peculiar momenta and T is the kinetic temperature. 
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For systems at equilibrium the thermodynamic temperature appearing in Eq. (2.3) should strictly 
only be calculated in the thermodynamic limit. However, for nonequilibrium systems, as mentioned 
above the application of the limit iV ^ oo is not straightforward. 

Still, for sufhciently large values of A^, Eq. (|2.2|) can be interpreted as: 



1 "i:' A,^(Fe) = + 0{N-^) = -^^^ + 0{N-^), (2.4) 

iv ^ — ' ukb ukb-l 

i—l 

where by 0{N~^) we mean that - at least at equilibrium - the finite size corrections can be bounded 
by a function of the form CN^^ with C of order 1. Although this condition on the constant C is 
not mathematically precisely defined, we discuss in section 5 numerical experiments, which will give an 
indication of the magnitude of the 0(N~^) corrections in the relationship between the KY-dimension 
and the viscosity. In what follows when we write 0{N^^), we will always intend it to carry the particular 
meaning given by the above discussion. 

Since we are interested here in the linear (Fg) regime one can use that LN{Fe) = Ln + 0{F^) 
where Xat = Ljv(O) / is the linear transport coefficient. We note that the left hand side of Eq. (^ 
is not restricted to small fields, due to its dynamical origin, and that the right hand side is related for 
small fields to the usual entropy production rate per unit volume of Irreversible Thermodynamics. £3 

We now introduce the Kaplan- Yorke (KY) dimension. If N^y is the largest integer for which 
Nky n 

^i,N{Fe) > 0, the KY-dimension, Dky.nM for a finite system with a discrete Lyapunov spectrum, 

defined by0§: 



1=1 
is 



n - AT a- ^i=l ^i,N{Fe) 

UkY.N - I>lKY + 7T 7T7TT (^-5) 

\ANKY+l-N{i'e)\ 



where we have not indicated the A^-dependence of Nky- 



3. SMALL PHASE SPACE REDUCTION 

In case the phase space dimension reduction is smaller than one, Eq. (^|^) can be reduced to the exact 
equation: 

V K^N{Fe) = X^^n,NiFe){2dN - / - DKYN{Fe)) = , (3.1) 



where the minimum Lyapunov exponent, \min,N{Fe) — ^2dN-f,N{Fe)- From Eqs. (3.1) and ( p.2| ) we 
can then trivially calculate the linear (i.e. the limiting zero field) transport coefficient in linear response 
theory as, 

,. {2dN - f - DKY,N{Fe))Xmax,N{Fe)nkBT 

= i™„ ■ (3-2) 

A similar relation has been obtained for the periodic Lorentz gas on the ba sis of periodic orbit theory. 
We remark that from the point of view of linear response theory, i.e. Eq. ( t3.2| ), a phase space dimension 
reduction smaller than one occurs for any iV, including N — > c». In that case one can let — > cxd in Eq. 



(3.2) and obtain a new exact relation for the linear transport coefficients, L = limjv^oo Ln, equivalent 
to the Green-Kubo formulae. The corresponding expression for the KY-dimension of the steady state 
attractor for sufficiently small fields, is given by: 

^Here and in what follows we assume that the transport coefficients are even in Fe. 
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DKY,N{Fe)^2dN-f- 



LF^N 



(3.3) 



To obtain the Eqs. (3.2) and (3.3), we have used that for systems which are symplectic at equihbrium, 

max,NiFe) — Knax,N +0{F^) = 



(i.e. ah Hamiltonian equihbrium systems), one can write for smaU F^: A. 



For systems which satisfy the Conjugate Pairing RuleEj (CPR) the sum of each conjugate pair, 
i, i* = 2dN — / — i + 1 of Lyapunov exponents is 



0(F2), where A 



max,N 



X. 



■rnax.N 



(0) and A, 



min.N • 



.A, 



■7nin,N 



(0). 



Az,Ar(i^e) + A,;.,Ar(Fe) 



-2aN{Fe)V 



Vz. 



(3.4) 



kB{2dN - /)' 

Wc note that in nonequihbrium systems, the Conjugate Pairing Rule is expected to hold only in systems 
that are thermostatted homogeneously. In some systems there is numerical evidence that the maximal 
expoiients satisfy the Conjugate Pairing Rule, while the other pairs do not (that is Eq. (3.4) is true for 
i=l).E3 These systems are said to satisfy the weak Conjugate Pairing Rule (WCPR). By combining the 
Eqs. (|3.1[) and (3.4), one obtains for sufhciently small fields: 



D 



KY.N 



(Fe) 



1 



A. 



■N 



(Fe) 



(dN - f/2) ^ A, 
(Isl) into Eq. K^) and using Eq, 



(3.5) 



Substituting Eq, 
limiting KY-dimension, for sufficiently small fields 



■rain, 

(3.3), one obtains another expression for the 



D 



KY,N 



{Fe 



{dN - f/2) 



= 3 



A 



min.N 



(Fe) 



A, 



max.N 



(Fe 



0{Ft) 



(3.6) 



In Eqs. ( [3.E| ) and ( ^.6| ), we have chosen to use the maximal Lyapunov exponents as the conjugate pair in 
Eq. ( |3.4[), a nd therefore these equations are valid provided WCPR is obeyed. A similar looking formula, 
as Eq. ( |3.6| ) with the first two terms on the right hand side only, has been quoted in p4[ . 

As mentioned before, all the results in this section hold under the hypothesis that the phase space 
dimension reduction is smaller than unity. For this to be true for any given small F^, N is constricted 
to be of 0{F^^), or equivalently, for any given large N, Fe to be of 0(A^~^/^), as can be seen from Eq. 
(3.3), so that Fe and N are coupled. |^ 

However, from a general physical point of view, we would like to have a theory for Dky.n, which 
holds uniformly in N, i.e. with Eq. ( |3.3| ) valid for every N with a small but fixed F^, so that N and Fe 
are independent variables. Now, it is trivial to generalize Eq. (B.l) to the case in which the dimensional 
reduction is greater than one. In fact one then obtains that: 



2dN-f 



J2 ^NiFe) = XNr,v + l.N{Fe){2dN - f - DkyMFc)) 



2dN-f 



(Ai,Ar(Fe) — XNKY + l,N{Fe)). 



(3.7) 



i = NKY+2 



If we assume, as is usually done, that for sufficiently large N, X2dN - f- j — X2dN-f + 0{N^^) for fixed j 
not varying with N; then for any fixed dimensional reduction, Eq. (3.7) simply becomes Eq. (3.1) with a 
correction of 0(1) in TV. Keeping the phase space dimension reduction fixed as N increases still imposes a 

^One might think that a phase space dimension reduction of one could hardly have any practical consequence in a 
macroscopic system whose phase space dimension is of the order of 10^^. As we will discuss in Section 6, such a very small 
phase space dimension reduction is expected to occur under physically realizable conditions. 
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condition on Fe of the form discussed above. To better control the errors when 2dN — / — Nky becomes 
large {0{N)), a more careful treatment of the Eq. (3.1) is needed in order to obtain an expression 
uniform in N for the phase space dimension reduction. 



4. LYAPUNOV SPECTRUM FOR VERY LARGE N AND 
LARGE PHASE SPACE REDUCTION 

For large N one can indeed reformulate the definition of Dky.n in a more analytical way. Consider 
thereto the stepwise continuous function of a continuous variable Q < x < 2dN — f: Xn{x, Fe) = Ai^jv(-Fe) 
for i — 1 < X < i. If one introduces the integral: 

fLN{x,Fe)= f ~XN{y,Fe)dy, (4.1) 







Eq. (2.5) can be rewritten in the form: 



ilNiDKYM{Fe),Fe) = 0. (4.2) 

Since DKY,N{Fe) as well as jj,N{x,Fe) are expected to be extensive quantities, i.e. they are proportional 
to N for large N , it is natural to define the quantity 

= 2dN-f (^-^^ 

where 5M{Fe) is the dimension per effective degree of freedom {2dN — /) in phase space. We now use 
again that, as N grows, the difference A^^at — \i+i^N is expected to go to zero as N^"^ . This suggests 
a possible rescaling of the variable x in \iq{x,Fe) to define the function AAr(a;,i^e) = Xn{xN, F^), 
i.e. XN{x,Fe) = Xi^N{Fe) for 2dN-f < X < 2dN-f "^^^ function Xi^]y{Fe) can then be expected to 
be well approximated by a continuous function oi the variable x when N is sufficiently large. Thus 



rewriting Eq. (4.1) as 

fiNix,Fe)^ I XN{y,Fe)dy (4.4) 



where < x < 1, Eq. (4.2) is equivalent to 



flN{SN{Fe),Fe)=0 (4.5) 

where 5is[{Ff,) and ^iN{x, F^) can be assumed to be intensive quantities. We can make this more precise by 
the following crucial assumption on the nature of the function Ajv(?/, Ff.), i.e. of the Lyapunov spectrum: 
Smoothness Hypothesis: If N is sufficiently large, one can write: 

XNix,Fe)^lix,Fe) + 0{N-^) (4.6) 
with l{x, Fe) a smooth fun ction^ of x and Fe 



Eq. (4^), as Eq. (p.4|), is to be interpreted that for sufficiently large N, \XNix,Fe) — l{x,Fe)\ < 
C'N^^ with C of order 1. Similarly, a s me ntioned already above, all the 0(7V~^) terms appearing 
in forthcoming equations (see e.g. Eq. ( 4.12| )) have to be interpreted in this way and the constants 



obtained then can be directly expressed in term of C . 



We now define, in analogy with Eqs. ([l.4|) and (4.5) 



is enough to assume that l{x, Fe) G in x and near x = 1. rMoreover we will require that l{x, Fe) is C'^ in Fe 
for every x. This may not be the case when a phase transition occurs. 113 
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Figure 1: Sketch of the integral of a continuous Lyapunov spectrum m{x,Fe) for a dissipative system, 
starting with the largest Lyapunov exponent at x — and ending with the smallest Lyapunov exponent 
at a; = 1. The dashed line is the tangent at x = 1. 



m{x,Fe) = / l{y,Fe)dy 
Jo 



and d{Fe) through the equation: 



m{d{F,),F,) = 
respectively. Clearly our Hypothesis implies that: 

6NiFe)^d{F,) + OiN-^) 



(4.7) 



(4.8) 



(4.9) 



where d{Fe) and m{Fe) are now intensive quantities. The function m(x,Fe) is sketched in Fig. 1 
especially near x — 1. 

From this figure one easily deduces that: 



d(Fe) = 1 



m(l,Fe) 



0(TO(l,i^e)') = 1 



1 



o(f4), 



(4.10) 



m'(l,Fe) ' V ' «^ / 1(1, Fe) 2dnkB 

where m'{Ff,) is the derivative of iin{x,Fe) with respect to x at x = 1. Here we used that m'{\,Fe) = 
Z(l, Fg) = Xmin{Fe) and that m(l, Fe)^ is the sum over all Lyapunov exponents, divided by {2dN — /), 
i.e. the phase space contraction (or entropy production) rate per (effective) degree of freedom in phase 
space.|^ Using Eq. (2^), Eq. ( 4.10 ) can be rewritten as: 



d[F,) = 1 



LF'i 



0{F^) 



(4.11) 



2dnkBT 

where the maximum and minimum Lyapunov exponents at equilibrium, (i.e. when F,, — 0) are Xmax — 
Z(0,0) and Xmin — '(1,0), respectively. Moreover we have used that Xmax — —Xmin and that l{\,Fi,) = 

?'2 



Z(1,0) + 0{F^). In terms of the extensive quantity DKY,N{Fe), Eq. ( 4.11 ) can be rewritten as 



F>KY.N{Fe) 



= 1 - 



LF? 



+ 0{Ft) + 0{N-'). 



(4.12) 



2dN - / X.,nax 2dnkBT 

where we kept the term f/{2dN) of 0{N-^ ) on the left hand side of Eq. ( |4.12|) to facihtate comparison 
in figure s 4 and 5 of section 5 for = 32 particles. We observe that Eq. ( 4.12 ) is formally very similar 
to Eq. (3^), except for the correction term 0{N~^) and the substitution of the asymptotic Xmax for 



^In a more analytic way Eq. (1.10) follows from the inverse function theorem together with our Smoothness Hypothesis 
and the fact that m(0, 0) = 0. 
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the finite N value Xmax^N- It is clear that the same generalization can be performed for the Eqs. 
and (3.6), if one notes that in the present context, the CPR, Eq. (|3.4|), becomes: 



l{x,F,) + l{l-x,F,) = - 



dnkn 



(4.13) 



From Eq. (4.11) it follows that in the linear regime the reduction in phase space dimension in large 
thermostatted dissipative systems is extensive. For those systems for which the Smoothne ss H vnothesis 
holds, this result is exact. The extensive nature of the reduction has been noted before. 



5. NUMERICAL TEST 



We tested our Smoothness Hypothesis as well as Eq. ( 4.12^ ) and equations derived from it assuming 
that the WCPR is valid, for a system of 32 WCA particleslH undergoing shear flow in two dimensions. 
Although_this system does not satisfy the CPR, it does appear to satisfy WCPR to within 0.7% when 
N — 32v3 (which was within the limits of numerical error adiieved). The equations of motion for 
particles in such a system are the so-called SLLOD equations, O 

cii =pi/m + ijyi, Pi ^Fi ~ ijpyi - api. (5.1) 

Here, at not too large Reynolds numbers, the momenta, p; are peculiar momenta, i is the unit vector 
in the x direction and Fi is the force exertedjm particle i by all the other particles, due to a Weeks- 
Chandler-Andersen pair interaction potentialO between the particles. The value of a is determined 
using Gauss' Principle of Least Constraint to keep the kinetic temperature fixed. H The SLLOD equa- 



tions of motion given by Eq.(5.1) then model Couette flow when the Reynolds number is sufhciently 
small so that laminar flow is stable. H f\ 

For this system the dissipative flux J is just the xy element of the pressure tensor Pxy] the transpoit 
coefficient L is the Newtonian shear viscosity 77; and the external field Fe is the shear rate 7 = du^ / dy 
where is the local flow velocity in the x-direction in the system. The calculations were carried out 
at a reduced kinetic temperature of unity and a reduced density n= N/V = 0.8. All physical quantities 



occurring in the Eq. (5.1) as well as the temperature and density are made dimensionless by reducing 
them with appropriate combinations of molecular quantities. In particular the reduction factor for 7 
amounts to about 1 ps~^ for Argon, (c.f. ||l^]) We note that for the SLLOD equatiQris_£or shear flow in 
2 dimensions using non autonomous Lees-Edwards periodic boundary conditions, O Ell / = 5 (due to 
the conservation of kinetic energy, momentum and position of the center of mass) . 

Figures 2 and 3 represent a direct test of our Smoothness Hypothesis. Figure 2 gives the Lyapunov 
spectrum for a system at a shear rate 7 = 0.05 and for a number of particles N = 8, 18, 32. As can 
be easily observed the Lyapunov exponents for iV=18 and 32 just fill the "open spaces" left by the 
exponents for A^=8 and 18, respectively. 

Figure 3 shows the behavior of Xmaxi^) and Xminil) as functions of the externally applied field 7. 
Although no numerical experiment can in general confirm a mathematical hypothesis, the numerical 
results seem to agree very well with our Smoothness Hypothesis. 

We now use the Lyapunov spectrum to compute Dky.n{i)- This is shown in figure 4a where the 
Dky.n{i) is plotted for N = 32, d — 2 and < 7 < 0.5. As can be easily seen, although from the 



definition Eq. (2.5) one would generally expect to see discontinuities in the first derivatives of this 
function for those values of 7 for which Nky change (c.f. |2^), the function appears very smooth for 
a rather large range of values of 7. This indicates that = 32 is already "big enough" to consider 

®The a appearing in the Eq. ( ^.l| ) is related to the phase space contraction rate A in Eq. ( |2.l[ ). It is given by the 
relation J^'^^^'i'^ A,; = A = -dNa + 0(1). 
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Figure 2: Discrete Lyapunov spectra for the SLLOD Eqs. ( p.l| ) for =8 (crosses), 18 (open circles) 
and 32 (small filled circles) and a shear rate 7 = 0.05, as a function of the scaled Lyapunov exponent 
pair index (i - 1)/{2N - 2), which runs from to 1 for 1 < i < 2A^ - 1. 
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Figure 3: Maximum and minimum Lyapunov exponents, Xmaxil) (crosses) and Xminil) (circles), re- 
spectively, as a function of 7. 
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2dN-f 




Figure 4: Comparison of various expressions for the KY-dimension per effective degree of freedom, 
Dky,n{'j)/ i'^dN — /) for sheared systems of = 32 particles in d = 2 plotted as a function of the 
shear rate 7 up to a st rain rate of a) 7 = 0.5 an d b) 7 = 0.2. Plotted are: D ky. n{i) I {"^dN — f) from 
its definition Eq. (2_^) (filled circle); from Eq. (5_^) ( plus sign); from Eq. (5_^) (crosses); and from 
Eq. (5.4) (open circle). The solid line is a fit to Eq. ( ^.5[ ) with a fourth order polynomial even in 7, 
and the dashed line identifies a phase space contraction of unity. We emphasize that in Eq. (5.2) the 



measured rj, deduced independently from the constitutive equation, has been used in the computation 
of Dky,n{i)- 



t he L yapunov spectrum as "effectively c ontin uous" .F[ This also permits us to check the validity of Eq. 
( 4.12 ), and the large N versions of Eqs. (3_^) and (|3^6|) , which in this case take the form: 



D 



KY.N 



(7) 



1 



riNT 



2dN - f 

Dky,n{i) 
dN - //2 

Dky.n{i) 
dN - f/2 



^rnax,N 



2dnkBT 

(7) 



■ N 



= 3 



..N 



(7) 
(7) 



A 



7nax,N 



(7) 



h 0{j^) + 0{N-^) 
f 0(7'*) + 0(A^-i) 
0{j^) + 0{N-'), 



(5.2) 
(5.3) 
(5.4) 



respectively. Figure 4a shows that these expressions can describe the system studied over the range of 
fields considered and that the correction terms are small. The deviations in the values obtained using 
the various expressions is within the nu meri cal error in the data. 

We observe here that Dky,n of Eq. (2^) is well approximated by a fourth order even polynomial in 
7 for fields up to 7 0.5 and that at 7 ^ 0.3 the terms of 0(7*) are just 5% of the terms of order 0(7^), 
while the phase space dimension reduction is clearly greater than unity. We note that the quadratic 
dependence of Dky,n on 7 implies a linear dependence of the shear stress on 7,123 so the linear regime 
for A'^ = 32 extends beyond strain rates where the phase space dimension reduction is smaller than one, 
which runs (for N = 32) till approximately 7 = 0.175. At 7 = 0.3, 77(7) is just 8% smaller than 7y(0). 

For completeness, in Figure 4b, we ex pan d Figure 4a in the regime where the phase space reduction is 
smaller than unity. In this regime, Eqs. (3^), (3_^) and (^^) are expected to be valid for the calculation 
of Dky,n and any deviations of the data from the value calculated using (2.5) are due to the limited 



''More precisely one can say that for TV = 32 the numerical errors involved in computing the Lyapunov exponents are 
already larger than the errors introduced by neglecting the 0{N^^) correction. This observation provides a more precise 
meaning of expressions like "sufficiently large TV" or "a constant C of order 1" . 



9 



D, 



0.997 
0.996 
0.995 
0.994 



2dN-f Q 

0.992 
0.991 
0.990 



2dN-f 



10 



20 



30 
N 



40 



50 



60 




Figure 5: Comparison of various expressions for the KY-dimension per effective degree of freedom, 
D ky.n{i) I {"^dN — f) for sheared systems of WCA particles in d = 2 at a strain rate of a) 0.15 and 
b) 0.5, plotted as a function of the number of particles, N. Plotted are: D ky,n {i)/ {'^dN — f) from 
(2.5) (filled circle); from Eq 

Eq, 



its definition Eq. 



|5^ ) (plus sign); from Eq. (5.3) (crosses); and from 
(3.3) can be used to approximate the system size at which as phase space 



Eq. (5.4) (open circle) 

dimension reduction of unity will occur for any given field. Using this relation, a phase space dimension 
reduction less than unity is expected to occur for system sizes smaller than that indicated by the dashed 
line. 



numerical precision of the independent calculations of the viscosity and the Lyapunov exponents, 0{F^) 
corrections, or to the assumption that the WCPR is obeyed. The numerical results indicate that the 
numerical error gives the most significant contribution to the deviations observed. The WCPR is found 
to be valid at least to within numerical error for the state points considered. This is consistent with 
previous work.O 

Using the fourth order fit to Dky.n, mentioned above, i.e. Dky,n{i) ~ 2dN — / + ^Dj^y + 
jDj^Y for the shear rate range [—0.2, +0.2], (c.f. figure 4b), one can calculate the zero strain rate 
or Newtonian shear viscosity from r/ — D N^^axnUsT /2N . Here Dj^y n ^^'^ ^ky n second 
and fourth derivatives of Dky.nIj), respectively, with respect to 7, taken at 7 = 0. This leads to a 
value of ?7 = 2.54 ± 0.07, which agrees very well within statistical uncertainties with the Newtonian 
viscosity directly measured in the simulation and calculated from the defining constitutive relation 
77(7) = —Pxy{l)h, viz. -q — 2.52 ± 0.05 (see the + points in figure 4b). Note that although the 
precision of Dky,n is high (less than 0.03% statistical error), the viscosity is related to the phase space 
contraction {2dN — f — Dky.n) and thus calculation of the viscosity from the phase space dimension 
reduction involves, at small fields, a very small difference between two large numbers, resulting in a larger 
relative error in the viscosity. Furthermore since the dimensional contraction is a quadratic function of 
the shear rate, it results in a more difficult calculation of the viscosity from the phase space dimension 
reduction than the constitutive equation which is linear in the shear rate at low strain rates. 

Figure 5 compares the values of Dky.n determined from Eqs. (5.2), (5.3) and ( ^.4| ) with those 
calculated from its definition in Eq. (p7^), as a function of N for fields of 7 = 0.15 and 7 — 0.5. 
In Figure 5a the results for the field of 7 = 0.15 are shown, and at this strain rate the dimensional 
contraction will be less than unity for systems size up to approximately N = 50. Therefore, for N ~ 50, 
Eq. (3^), which contains no 0{N~^) corrections, will be valid, as will Eqs. ( |3.5| ) and ( |3.6| ) if the WCPR 
is obeyed. All the|_numerical results are consistent with the theory and with the assumption that the 
WCPR is obeyed. Il3 At 7 = 0.15, the values determined using the various methods have a maximum 
difference of 0.2% and are within the numerical errors at each particle number. In Figure 5b, the results 
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for a strain rate of 7 = 0.5 are shown. In this case, a dimensional contraction of less than unity is only 
obtained for iV ~ 5. Again, the deviations between the results calculated using Eqs. (5.2)- (5.4) and 
the definition of Dky,n given by Eq. ( |2.5| ), are small (at most 1%), and within the limits of error for 
all particle numbers considered. This confirms that the coefficients of the 0{N^^) and O(F^) terms are 
small for this system. 



6. CONCLUSIONS 

We mention here a few implications of the results presented in this paper. 

1. The extensivity of the phase space reduction for large N and small fields is here, to the best of 
our knowledge, demonstrated for the first time, on the basis of the Smoothness Hypothesis of the 
Lyapunov spectrum and the extensivity of the total entropy production. 



The relationships given by Eqs. (3.2), (3.3) and (5.2) also apply to systems where not all particles 



are thermostatted. That is, they can be applied to systems where the Gaussian thermostat operates 
on selected particles, say those in the boundaries, while the jmmaining particles evolve under 
Newtonian dynamics, supplemented by a dissipative field. O c3 We note that the Eqs. (|3.5[), 



(|3.6D, (5.3) and (5.4) can only be assumed to apply to homogeneously thermostatted systems in 



general, since only for such systems can the WCPR be expected to hold. 

A simple calculation shows that for a typical case as that of one mole of Argon at its triple 
point, sheared at the rate of 1 Hz, the difference of the Kaplan- Yorke dimension and the phase 
space dimension (0(10^'^)) is tiny, namely ^ 3. This follows from Eq.(3.3), which shows that the 
dimension loss, when measured in moles, is equal to the product of the total entropy production 
rate of the system and the reciprocal of the largest Lyapunov exponent. Since the largest Lyapunov 
exponent is controlled by the most unstable atomic processes, it is always very small ^ 1 ps^^ 
whether for atomic, molecular and even polymeric systems. 

We note that this smallness of the phase space dimension reduction in irreversible processes near 
equilibrium could well be the reason that linear Irreversible Thermodynamics provides such a good 
description of nonequilibrium systems close to equilibrium. This is because the thermodynamic 
properties are insensitive to the high order distribution functions - including the full N-particle 
distribution function of the entire system - since they are determined by a few low order distribution 
functions, which "do not know" that the dimension of the steady state attractor is only a few 
dimensions smaller than the ~ 10^^ of the phase space of the system. 



For the system studied, the 0{N ^) corrections to Dky,n / {"^-dN — /) that appear in Eqs. ( ^.2| 



( p.3|) and ( |5.4|) , due to the smoothness hypothesis are small (see Figure 5) and less than 1.0% even 



for system sizes of = 6. 
In conclusion, although the new relations involving Dky,n are simple consequences of the Eq. 



(4.10), they could nevertheless be useful for applications. In particular the equations give a simple 
expression for the Kaplan- Yorke dimension of the attractor of a class of many particle systems 
close to equilibrium i.e. in the regime of linear dissipation near equilibrium. 
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